Introduction to the Bethe Ansatz III 
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Having introduced the magnon in part I and the spinon in part II as the relevant quasi-particles 
for the interpretation of the spectrum of low-lying excitations in the one-dimensional (ID) s=l/2 
Heisenberg ferromagnet and antiferromagnet, respectively, we now study the low-lying excitations of 
the Heisenberg antiferromagnet in a magnetic field and interpret these collective states as composites 
of quasi-particles from a different species . We employ the Bethe ansatz to calculate matrix elements 
and show how the results of such a calculation can be used to predict lineshapes for neutron scattering 
experiments on quasi-lD antiferromagnetic compounds. The paper is designed as a tutorial for 
beginning graduate students. It includes 11 problems for further study. 
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I. INTRODUCTION 

In most areas of condensed-matter research, a few 
model systems outshine all others by their prototypical 
significance, which promises to encapsulate the essence 
of a physical phenomenon with little interference of col- 
lateral degrees of freedom. In cooperative magnetism, 
the Ising model, the Heisenberg model, and the Hubbard 
model, for example, are theoretical many-body systems 
that have kept generations of researchers fascinated. 

The prominence of such models and the abundance of 
theoretical predictions for many of their physical prop- 
erties have at times led to a reversal of the traditional 
relationship between theory and experiment. Instead of 
theorists proposing and solving models for the purpose 
of explaining measurements performed on specific ma- 
terials, we see magneto-chemists and condensed-matter 
experimentalists at work searching for materials which 
are physical realizations of the prototypical models. 

The focus here is on the one-dimensional (ID) s = 
i Heisenberg antiferromagnet in a magnetic field. The 
Hamiltonian for a cyclic chain with N sites reads 



N 

H = [JS„ • S„ + i — hS^] 

71=1 



(1) 



This model is amenable to exact analysis via Bethe 
ansatz and displays dynamical properties of intriguing 
complexity. The magnetic field h is a continuous parame- 
ter which has a strong impact on most physical properties 
and which is directly controllable in experiments. The 
high degrees of computational and experimental control 
is what makes this system so attractive to researchers. 

Physical realizations of Heisenberg antiferromagnetic 
chains have been known for many years in the form of 
3D crystalline compounds with quasi-lD exchange cou- 
pling between magnetic ions. The desired properties of 
the best candidate material include the following: The 
coupling of the effective electron spins must be highly 
isotropic in spin space and overwhelmingly predominant 
between nearest-neighbor magnetic ions along one crys- 
tallographic axis. The intra-chain coupling must not 



be too weak or else it will be hard to study the low- 
temperature properties, which are of particular interest. 
It must not be too strong either or else it will be hard 
to reach a magnetic field that makes the Zeeman energy 
hS^ comparable to the exchange energy JS n ■ S n +i. 

One compound that fits the bill particularly well is 
copper pyrazine dinitrate [Cu(C4H4N2)(N03)2], a mate- 
rial that was only recently synthesized in single crystals of 
sizes and shapes suitable for magnetic inelastic neutron 
scattering. EJ The most detailed dynamical experimental 
results to date were made on KCuFap in which the intra- 
chain coupling is considerably stronger, which makes it 
easier to study low-temperature effects but harder to 
study magnetic-field effects. 

What happens in a magnetic neutron scattering ex- 
periment? A beam of monochromatic neutrons, ideally a 
plane wave with well defined momentum and energy, is 
scattered inelastically off an array of spin chains via mag- 
netic dipolar interaction between the exchange coupled 
electron spins and the neutron spin into a superposition 
of waves with a range of momenta and energies. Each in- 
elastic scattering event causes a transition between two 
eigenstates of the spin chain. The difference in energy 
and wave number of the two eigenstates involved must 
be matched by the energy and momentum transfer of 
the scattered neutron. 

If the measurement is performed at very low tempera- 
ture, it is reasonable to assume that the observable scat- 
tering events predominantly involve transitions from the 
ground state \G) to excitations |A) with energies and 
wave numbers within a window preset by the experimen- 
tal setup. The experiment thus enables us to have a 
direct look at the excitation spectrum of the spin chain. 
It probes the spin fluctuations as described by the fluc- 
tuation operator 



N 



N- 1/2 J2 el9nS n> 



fi = x,y,z 



(2) 



for wave numbers q = 2irl/N, 1 = 1,.., ,N. Under ideal- 
ized circumstances, the inelastic neutron scattering cross 
section is proportional to the dynamic spin structure fac- 
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tor at zero temperatureiaa 

S w (g, w) = 2tt^ I (G\S%\\) \ 2 S (w-wx). (3) 

A 

Each scattering event with energy transfer huj\ = 
E\ — Eq and momentum transfer hq = h(k\ — ka) 
along the chain induces a transition from \G) to |A) and 
contributes a spectral line of intensity 2tt\ (G\Sjf\\) | 2 to 
Sfj,ij,(q, u>). For a macroscopic system, the spectral lines 
are arranged in a variety of patterns in (q,w)-space, in- 
cluding branches, continua, and more complicated struc- 
tures. The functions S^(q,w) thus provide a palette of 
information about the T = dynamics of the spin chain. 

Unlike in a classical dynamical system, where all mo- 
tion grinds to a halt at T = 0, the quantum wheels keep 
turning - a phenomenon known as zero point motion. 
For the Heisenberg model (H) in the ground state \G), 
the quantum fluctuations of any observable O of interest 
can be described by a time-dependent correlation func- 
tion of the form 

(0(t)O f ) ee (G\e tHt/h Oe~ lHt/h O^\G). (4) 
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FIG. 1: Energy versus wave number of all 2-spinon scattering 
states with S T — St — 1 at q > for TV = 64 (o) in com- 
parison with the corresponding free 2-spinon superpositions 
(+). The large red circles are 2-spinon data for N = 16. The 
inset shows the energy-momentum relation ([]) of the spinon 
quasi-particle. 



If O is conserved ([0,H] = 0), then (0(t)0^) = 
(G\00'\G) is a constant. Different dynamical vari- 
ables O yield different correlation functions for the same 
ground state \G). The dynamical variable seen by neu- 
trons during the magnetic scattering is S%. Fourier trans- 
forming (|j) with O — S^ yields the dynamic spin struc- 
ture factor (||) (Problem 1). 

Light scattering, electron scattering, photoemission, 
and nuclear magnetic resonance are some other exper- 
imental techniques used to investigate the dynamics of 
(|l|). Different measuring techniques view the same exci- 
tation spectrum through lenses of different color, i.e. by 
transition rates specific to particular fluctuation opera- 
tors. Hence each experimental probe filters out a specific 
aspect of the zero point motion by viewing a particular 
dynamical variable of one and the same system. 

The goal set for this column is to teach the reader how 
to make detailed predictions based on the Bethe ansatz 
solution of ([!]) for the inelastic neutron scattering cross 
section as can be measured on quasi-lD antiferromag- 
netic materials. Many of the computational and analytic 
tools tkat are needed for this purpose were introduced in 
parts IB and HQ of this series. The calculation of matrix 
elements from Bethe wave functions will be introduced 
here along the way. 



II. QUASI-PARTICLES AT h = 0: SPINONS 

In part II we used the Bethe ansatz to describe the 
ground state \A) of (0) in zero field. We saw that this 
singlet (St = 0) state can be interpreted as the physi- 
cal vacuum of a. particular species of particles with spin 
i: the spinons.EI To distinguish them from the electrons, 



protons, and neutrons - the constituent elementary par- 
ticles of all materials - we use the name quasi-particle 
for the spinons as is custom. We identified a set of low- 
lying excitated states containing two spinons with spins 
up. The spin S T = St = 1 of a (stationary) 2-spinon 
state is shared by all electrons of the system in what is 
called a collective excitation. 

Here we revisit these 2-spinon excitations with our eyes 
focused on the quasi-particles. The red circles in Fig. |l| 
represent the 2-spinon states for N = 16 (see also Fig. 4 
of part II). We know that for N — > oo they close-, up to 
form a continuum in (q, w)-space with boundaries&E 



tu(q) = 



■ q 

sm- 



(5) 



represented by the solid lines in Fig. [l]. In every eigen- 
state of this set, the spinons can be thought of as two 
localized perturbations of the spinon vacuum \A) mov- 
ing around the chain with momenta pi , P2 and energies 
e sp(Pi), £sp(p2)- The spinon energy-momentum relation, 



e s P (p) = -z Jsin P' 0<p<7r, 



(6) 



is shown in the inset to Fig. |l|. Periodically, the two quasi- 
particles scatter off each other, hence the name 2-spinon 
scattering state. 

The energies and wave numbers of the 2-spinon eigen- 
states for N = 64 obtained via Bethe ansatz are repre- 
sented by circles (o) in Fig. |l], whereas the corresponding 
(fictitious) free 2-spinon superpositions with wave num- 
bers q = pi + p 2 and energy u) = e sp (pi) + £ S p(p2) are 
shown as (+) symbols. The upward displacement of the 
former relative to the latter is a measure of the positive 
spinon interaction energy caused by a predominently re- 
pulsive force acting between the two quasi-particles. 
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FIG. 2: Scaled spinon interaction energy e ap (p l ,p2) in the 
2-spinon states for N = 64 versus spinon momenta pi , p2 (3D 
plot) and versus relative spinon momentum \p\ — p%\ (inset). 

With increasing N, the scattering events between the 
two spinons in a 2-spinon state become scarcer at a 
rate inversely proportional to the distance covered by 
the quasi-particles between collisions. Consequently, the 
spinon interaction energy is expected to diminish oc N^ 1 
(Problem 2). We conjecture that, asymptotically for 
large N, the spinon interaction energy in the 2-spinon 
scattering states, as defined by the expression 

AE^iq) ee E^(q) - e sp (p x ) ~ (?) 
with q = pi + p2 has the form 

AE^ p \q)=e sp ( Pl ,p 2 )/N. (8) 

The quantity e sp {pi,P2) can then be attributed to a single 
2-body collision between spinon quasi-particles traveling 
with momenta p\,P2- 

The N — 64 data for e sp (p\,p2) shown in Fig. || indi- 
cate that the interaction energy depends smoothly on the 
spinon momenta. It is observed to be smallest when both 
spinons have equal momenta. These states are near the 
upper boundary of the 2-spinon continuum. For fixed 
momentum p\ > P2, the interaction energy is a mono- 
tonically decreasing function of p2- Hence the largest 
interaction energy of states at fixed wave number q is re- 
alized when the relative spinon momentum \p% — P2\ is a 
maximum. These states are near the lower boundary of 
the 2-spinon continuum. 

At fixed \p\ — p%\ , the spinon interaction energy de- 
pends only weakly on the wave numbers of the 2-spinon 
state. Glancing at the data in Fig. g along the lines 
Pi + Vi =const. (direction of arrow), makes them col- 
lapse into a narrow band in the plane spanned by e sp 
and \p\ — p2\ as shown in the inset to Fig. ^. 

In Fig. 3 of part I we had illustrated the magnon in- 
teraction energy for 2-magnon scattering states. There 



the smallest interaction energy was realized when one of 
two magnons had infinite wavelength {j>2 =0). It then 
acted like a slightly rotated magnon vacuum in which the 
other magnon could move freely. A spinon, by contrast, 
becomes the strongest scatterer to another spinon when 
it has zero momentum. 

The vanishing spinon interaction energy for N — > oo 
does not make the calculation of transition rates for the 
dynamic structure factor (||) a simple task. It was only 
recently that the exact 2-spinon part o£-S zz (q, lo) at h = 
was evaluated for an infinite chain.£3 The techniques 
used in that calculation, which involved generating the 2- 
spinon states from the spinon vacuum by means of spinon 
creation operators and expressing the spin fluctuation op- 
erator (|^) in terms of spinon creation operators, are not 
readily generalizable to h ^ (Problem 3a). 

In nonzero field, a different approach is required. The 
ground state must be reinterpreted as the physical vac- 
uum for a different species of quasi-particles. The spec- 
tral weight of S zz (q,ui) will be dominated by scattering 
states of few quasi-particles from the new species, and the 
associated transition rates will be calculated via Bethe 
ansatz. 



III. QUASI-PARTICLES AT h 0: PSINONS 

Increasing the magnetic field from ft, = to the satu- 
ration value h = hs — 2J, leaves all eigenvectors of (|l|) 
unaltered but shifts their energies at different rates. The 
ground state |G) changes in a sequence of level cross- 
ings, and the magnetization increases from M z = to 
M z = N/2 in units of one as discussed in part II. Every 
level crossing of that sequence adds two spinons to \G). 

At h the number of spinons moving through the 
chain and scattering off each other is of O(N). Scattering 
events are frequent. The spinon interaction energy does 
not vanish in the limit N — ► oo, because the density 
of spinons remains finite. Hence the spectrum of low- 
lying excitations cannot be reconstructed from the spinon 
energy-momentum relation (||) . 

If we start from the ground state \F) = | TT " " T) a t 
h — hs instead and decrease the the magnetic field grad- 
ually, we are also facing the problem of a runaway quasi- 
particle population. The state \F) is the vacuum of 
magnons, a species of quasi-particles with spin 1. Every 
reverse level crossing of the sequence considered previ- 
ously adds one magnon as it removes two spinons. In 
the 2-magnon scattering states, which we had studied 
in part I, the role of individual magnons was just as 
clearly recognizable as the role of individual spinons in 
the 2-spinon scattering states depicted in Fig. [|. But 
in \G) at < h < hs and in all low- lying excitations 
from that state, the number of magnons is macroscopic 
(oc N). The magnon interaction energy remains signifi- 
cant, which makes it impossible to infer spectral proper- 
ties from the magnon energy-momentum relation (15). 

The ground-state wave function \G) at < h < hs has 
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FIG. 3: Psinon vacuum \G) for a chain of N = 8 spins at 
magnetization M z = 0, 1, . . . ,4. The position of the magnons 
(•) are determined by the set (^[) of Ii's and the positions of 
the spinons (Q) by the vacancies across the full range of the 
ii's allowed by (0). 

spin quantum numbers S T = S% = M z (0 < M z < N/2) 
and is specified by the following set of r = N/2 — M z 
Bethe quantum numbersO 



{Ii}a - \ -j + -7T + o' • • • > 7T~ * 



M, 
~2 



N 

T 



M z 
~2 



(9) 



Henceforth we treat \G), which can be interpreted as a 
state containing N/2 — M z magnons or as a state con- 
taining 2M Z spinons, as a new physical vacuum and 
describe the dynamically relevant excitation spectrum 
from this reference state by modifying the uniform array 
(0) of Bethe quantum numbers systematically. For this 
purpose we consider the class K r of eigenstates whose 
Bethe quantum numbers comprise, for < r < N/2 and 
< to < N/2 — r, all configurations 



m < 1% < h < ■ ■ ■ < I r < 



(10) 



We recall from part II that the i,'s must be integer valued 
if r is odd and half-integer valued if r is even. Every class- 
ier eigenstate is represented by a real solution {zi} of the 
Bethe ansatz equation (115), 



N0(z t ) = 2nI l +Yj[(2 



Zj )/2], i=l,. 



(11) 



with 4>{z) = 2arctanz. These solutions can be obtained 
iteratively via (119) for fairly large systems. Every class- 
ier state at fixed integer quantum number m (0 < to < 
M z ) can be regarded as a scattering state of m pairs of 
spinon-like particles, which we have named psinons. 

At M z = the psinon vacuum coincides with the 
spinon vacuum, both containing N/2 magnons. At sat- 
uration [M z = N/2), the psinon vacuum coincides with 
the magnon vacuum, a state containing N spinons. The 
transformation of the psinon vacuum between the spinon 
vacuum and the magnon vacuum is illustrated in Fig. |^. 
It displays the configuration of Bethe quantum numbers 
for \G) in a system with N = 8 and all values of M z 
realized between h = and h = hs- 

The top row (spinon vacuum) corresponds to the top 
row in Fig. 3 of part II (albeit for different N) . The sec- 
ond row is a particular state of the 2-spinon triplet set 
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FIG. 4: Psinon vacuum (m = 0) for N = 8, M z = 2 and 
complete set of 2m-psinon excitations with m = 1, 2 and wave 
numbers q = k — he (in units of 2-k/N). The positions of the 
magnons (small circles) are determined by the Ii's and the 
positions of the spinons (large circles) by the U vacancies. A 
subset of the spinons are the psinons (green circles). 



discussed in part II. The third row represents the psinon 
vacuum at half the saturation magnetization, M z /N = i, 
the case we shall focus on in this paper. Here the psi- 
non vacuum contains twice as many spinons as it con- 
tains magnons. The psinon vacuum in the fourth row 
corresponds to the highest-energy state of the 1-magnon 
branch described in part I. Trading the last magnon for a 
pair of spinons saturates M z and makes the psinon vac- 
uum equal to the magnon vacuum (fifth row). 

In the psinon vacuum \G), the only class- K r state with 
to = 0, the magnons form a single array flanked by 
two arrays of spinons. Relaxing the constraint in ( |To| ) 
to to — 1 yields the 2-psinon excitations. Generically, 
magnons now break into three clusters separated by the 
two innermost spinons, which now assume the role of psi- 
nons. The remaining 2M Z — 2 spinons stay sidelined. 

In the 4-psinon states (to = 2), two additional spinons 
have been mobilized into psinons. By this prescription, 
we can systematically generate sets of 2m-psinon exci- 
tations for < to < M z . From the psinon vacuum for 
N = 8, M z = N/4 = 2 depicted in the the third row of 
Fig. H we can thus generate the five 2-psinon states and 
the nine 4-psinon states identified in Fig. Q 

Figure ^| shows energy versus wave number of all 2- 
psinon states at M z /N = | for system sizes N = 16 
(large red circles) and N — 64 (small circles). In the 
limit N — ► oo, the 2-psinon states form a continuum in 
(q,oj)-space with boundaries shown as blue lines. The 
2-psinon spectrum is confined to \q\ < q s , where 



tt(1 - 2M Z /N). 



(12) 
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FIG. 5: Spectrum in (q, tj)-space of all 2-psinon scattering 
states at M z /N — j for TV = 16 (large red circles) and TV — 64 
(small circles), the latter in comparison with the correspond- 
ing (fictitious) free 2-psinon states (+). The spectral range 
of the 2-psinon continuum for TV — > oo (blue lines) are in- 
ferred from data for TV = 2048. The inset shows the psinon 
energy-momentum relation. 

From the 2-psinon continuum boundaries we infer the 
psinon energy-momentum relation (p) shown in the in- 
set, by the requirement that the wave number and the en- 
ergy of every 2-psinon state for TV — > oo can be accounted 
for by q = pi + p 2 and oj = e^(p{) + e^(p 2 ), respectively. 
The two arcs that make up the lower 2-psinon continuum 
boundary are then given by e^(q ± 7r/4). 

For finite TV, the scattering of the two psinons in the 
2-psinon states produces an interaction energy. In Fig. ||, 
this energy can again be measured by comparing the 
positions of the TV = 64 scattering states (o) relative 
to the positions of the corresponding (fictitious) free 2- 
psinon states (+). The evidence from a comparison of 
the TV = 64 data in Figs. || and ^ is that the psinons at 
M z /N = I interact more strongly than the spinons at 
M z = in chains of the same length. The interaction 
energy between two psinon quasi-particles in a 2-psinon 
scattering state can be investigated by the method used 
m Sec. for spinon (Problem 4a). It again varies oc 1/TV. 

Mobilizing two additional spinons into psinons makes 
the 4-psinon scattering states much more numerous than 
the 2-psinon states. Their number grows oc TV 4 . In 
Fig. U we have plotted all 4-psinon states for TV = 16 
and TV = 64 as well as the 4-spinon spectral range. The 
4-psinon spectral threshold has the same shape as the 
2-psinon lower boundary but extended periodically over 
the entire Brillouin zone. The upper 4-psinon boundary 
is obtained from the upper 2-psinon boundary by a scale 
transformation (q — > 2q,to — > 2u>). 

In both sets of states, the energy correction due to 
the psinon interaction is a 1/TV effect. It fades away as 
the scattering events become less and less frequent in 



FIG. 6: Spectrum in (q, o;)-space of all 4-psinon excitations 
at M z = TV/4 for TV = 16 (red squares) and TV = 64 (dots). 
Also shown is the range of the 4-psinon states for TV — > oo 
(blue lines) inferred from TV = 2048 data. 

a system of increasing size. A comparison of finite-TV 
data in Figs. [| and [6] shows that the finite-size energy 
correction caused by the psinon interaction is stronger 
in the 4-psinon states than in the 2-psinon states. The 
reason is the different rates at which scattering events 
between psinons occur. 

At M z = the 2-spinon excitations were found to 
dominate the spectral. weight in the dynamic spin struc- 
ture factor S zz (q, w)£3 Our task here is to determine how 
the spectral weight of S zz (q,u>) at M z ^ is distributed 
among the 2m-psinon excitations. To accomplish it, we 
must calculate matrix elements. 



IV. MATRIX ELEMENTS 

The Bethe ansatz has rarely been used for the purpose 
of calculating matrix elements. Most attempts at taking 
this approach have been deterred by the need of evaluat- 
ing the sum V £ S r over the r! magnon permutations in 
the coefficients (128) of the Bethe eigenvectors (127): 

\ip) = a(ni,... ,n r )\ni,... ,n r ), (13) 

i<m<...<n r <iV 

a(m, ■ ■ . , n r ) = ^2 ex P * ^2 kv 3 n J + o ^ 0viVj ' 
ves r \ 3=1 i<j J 

The magnon momenta {k{\ and the phase angles {Oij} 
are related via (114) and (115) to the solutions {zi} of the 
Bethe ansatz equations (O). 

In the calculation of a single matrix element, a sum 
V G S r is evaluated many times, once for every coefficient 
a(ni,... ,n r ) of the two eigenvectors involved. Under 
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these circumstances, it is imperative that the algorithm 
has rapid access to a table of permutations. 

Table | describes a powerful recursive algorithm that 
generates all permutations of the numbers 1,2,... , r.L3 
These permutations are stored in the array ip of dimen- 
sionality r! x r. The storage of this array in the RAM 
requires 645kB for r = 8, 6.53MB for r = 9, and 72.6MB 
for r = 10. 

The computational effort of calculating any Bethe 
ansatz eigenvector ( |l3| ) can be reduced considerably if 
we use the translational symmetry, Tl^) = e tk \ip) (see 
part I). It is guaranteed by the relation 



a(ni + I, . . 



+ l) = e lkl a{n 



(14) 



between sets of coefficients pertaining to basis vectors 
that transform into each other under translation. Here 
the integers rij + 1 have to be used mod(A). Translation- 
ally invariant basis vectors have the form 



li;fc) = 



where = T z |j) 



1 



E e " fc i^ 



"0 l=Q 



I. 



,00 



(15) 



I) and 1 < 



N/dj < N is an integer. The wave numbers k realized in 



the set ( |15[ ) are multiples mod(27r) of 2-ir/dj 

M 
i 1 , . . 



The set of basis vectors \j)o — \n 
1, . . . , d, are the generators of the translationally invari- 
ant basis. The set of distinct vectors \j; k) for fixed k 
is labeled j € Jk Q {1, • • ■ , d}. The rotationally invari- 
ant subspace for fixed N/2 — r, which has dimensionality 



D = N\/[r\(N — r)!], splits into N translationally invari- 
ant subspaces of dimensionality Dk, one for each wave 
number k = 2im/N, n = 0, . . . , N — 1. We have 

d 

D = J2 d i= E D *> D * = E ■ ( 16 ) 

j=i o<fc<27r jej k 

The Bethe eigenvector ( p^ ) expanded in this basis can 
thus be written in the form 



lv>> = E a ;E ei "^ = E 

j£Jk 1=0 3 = 1 



a E eilk \fil> 

1=0 



(17) 



where the a,j = a(n^ , . . . ,ni. J '), the Bethe coefficients 
of the generator basis vectors |j)o, are the only ones that 
must be evaluated. The last expression of (0) holds 
because the Bethe coefficients aj of all generators \j)o 
which do not occur in the set Jk are zero (Problem 5). 

We calculate transition rates for the dynamic structure 
factor (||) in the form 



l(G|^|A)| 2 = 



\\u?\w 



(18) 



where |^q),|^x) are the (non-normalized) Bethe eigen- 
vectors of the ground state and of one of the excited 
states from classes (i)-(vi), respectively. The norms are 
evaluated as follows: 



E d ,i 



(19) 



TABLE I: C++ code of a recursive algorithm that generates 
the r! permutations of the numbers 1,2,... , r and writes them 
into an array of dimensionality r! xr. Each number is stored as 
a chaq or short , which is of size 2Bytes on an Intel processor 
(RefB). 



void main() 
{ 

perm(0,r) ; 

} 

void perm(int k, int r) 
-C 

static long id=-l; 
static long rf=l; 
ip[0] [k]=++id; 
if (id==r){ 

for (int l=0;Kr;l++) ip [rf ] [1] =ip [0] [1+1] ; 
rf++; 

} 

for(int t=l;t<=r;t++) if (ip [0] [t] ==0) perm(t.r); 
id—; 

ip[0] [k]=0; 

} 



The matrix element (i/)o\Sjf\ip\) is nonzero only if q = 
k\ — fco + 27rZ. For the fluctuation operator S* it can be 
evaluated in the form (Problem 6): 



1 



d N 

E4 0) 4 A) E 



dj-l 

"E 1 

1=0 



Ah. 



(j\S z n \j) 



(20) 



The non- vanishing matrix elements (■)/>o|<S'^ : |t/'a) needed 
for S xx (q, must also satisfy q — k\ — ko + 27rZ and can 
be reduced to somewhat more complicated expressions 
involving elements i (jo\S^\j\}i x between basis vectors 
from different subspaces. 

What are the memory requirements for the calcula- 
tion of one such matrix element? Three arrays are 
needed for the basis vectors: long b [d] , short nd [d] 
and, ndr [d] [r] . The array b [d] holds the bit pattern 
of \j)o, nd[d] holds the numbers dj, and in ndr[d] [r] 
we store the r numbers rv{, ... ,nl belonging to \j)o, i.e. 
b [ j ] . The arrays nd , ndr are not really necessary, since 
we can always calculate those numbers from b [ j ] , but 
they help reduce the CPU-time. Finally, the two arrays 
complex psi[d] hold the coefficients aj. Again, storing 
these numbers for multiple use cuts down on CPU time. 



7 



TABLE II: Memory (in kB) required for the calculation of 
matrix elements via Bethe ansatz on an Intel processor. 
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In Table we have summarized the memory require- 
ments for several applications to the Heisenberg model 
(|l|). With the exception of the case of N — 20, r — 10, we 
can compute the transition rates on a simple Pentium-PC 
with at least 32MB of memory. It is also worth mention- 
ing that for the case N = 16, r = 8 no more than 1MB 
of memory is needed. 



V. PSINONS AND ANTIPSINONS 

Which 2m-psinon excitations have the largest spec- 
tral weight in S zz (q,uj)7 We begin with a chain of 
N = 16 spins at magnetization M z /N = j, and explore 
the transition rates between the ground state |G) with 
{Ii\a = ^{ — 3, — 1, +1, +3}, and all 2m-psinon excita- 
tions for ui — 0,1, 2, 3, 4. 

First we calculate (G\Sq\G), which probes, for q = 0, 
the ground-state magnetization induced by the magnetic 
field (Problem 7). Next we investigate the 2-psinon 
states. The Ii configurations are shown in Fig. 0. The 
first row represents the psinon vacuum with its four 
magnons sandwiched by two sets of four spinons. The two 
innermost spinons (marked green) become psinons when 
at least one of them is moved to another position. In the 
rows underneath, the psinons are moved systematically 
across the array of magnons while the remaining spinons 
stay frozen in place. These eight configurations describe 
all 2-psinon states with wave number < q < tt/2. 

The lis of the states shown in Fig. [j] and their wave 
numbers as inferred from (II8b) are listed in Table |l|. 
Solving the Bethe ansatz equations ( pd] ) yields the in- 
gredients needed to evaluate the energies via (II8a) and 
the transition rates via (|20|), which are also listed in Ta- 



bic [II . We observe that almost all of the spectral weight 
is concentrated in the lowest excitation for any given 
wave number (lowest red circles in Fig. |^). In a macro- 
scopic system, these states form the lower boundary of 
the 2-psinon continuum. 

When we calculate the 4-spinon transition rates 
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FIG. 7: Psinon vacuum |G) for N = 16, M z = 4 and 2-psinon 
states with q > 0. The Ii values are marked by the positions 
of the magnons (small circles). The spinons (large circles) 
mark Zi-vacancies. A subset of the spinons are called psinons 
(green circles). 



TABLE III: Ground state and 2-psinon excitations for N = 
16, M z = 4, and wave numbers q = k — kc > (in 
units of 2-rr/N). The ground state has ka — and Eg = 
-11.5121346862. 
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| | A) | 2 , we find that most of their spectral weight is 
again carried by a single branch of excitations. The dy- 
namically dominant 4-psinon states for 7Y = 16 are the 
eight lowest red squares in Fig. ^. For large N they form 
a branch adjacent to the 2-psinon spectral threshold. An 
investigation of the 2m-psinon states for m = 3, 4, . . . 
shows that there exists one dynamically dominant branch 
of 2m-psinon excitations for < m < M z . All other 2m- 
psinon excitations have transition rates that are smaller 
by at least two orders of magnitude at q < tt/2 and still 
by more than one magnitude at q > tt/2. 

Figure |s| shows the Ii configurations of the four dynam- 
ically dominant 2™-psinon branches for N = 16. Each 
branch (at q > 0) consists of A^/2 — M z — 4 states. An 
interesting pattern emerges, which is indicative of the na- 
ture of the relevant quasi-particles in these dynamically 
dominant collective excitations. The two relevant quasi- 
particles are highlighted by green circles. 

We identify one of the two quasi-particles as a psinon 
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FIG. 8: Psinon vacuum \G) for N = 16, M z = 4 and set of 
V'V'* states with < q < it. The 7i are given by the positions 
of the magnons (small circles) in each row. The spinons (large 
circles) correspond to /^-vacancies. The psinon (if)) and the 
antipsinon (ip*) are marked by a large and a small grey circle, 
respectively. 



(large green circle) as before and the other one as a new 
quasi-particle (small green circle). The latter is repre- 
sented by a hole in what was one of two spinon arrays of 
the psinon vacuum. Instead of focusing on the cascade 
of psinons (mobile spinons) which this hole has knocked 
out of the psinon vacuum, we focus on the hole itself, 
which has properties commonly attributed to antiparti- 
cles. The psinon (ip) and the antipsinon (jp*) exist in dis- 
junct parts of the psinon vacuum, namely in the magnon 
and spinon arrays, respectively. When they cross paths 
at the border of the two arrays, they undergo a mutual 
annihilation, represented by the step from the second row 
to the top row in Fig. |^. 

The large red circles in Fig. ^| represent all ipip* states 
at q > for N — 16 as specified in Fig. ||. There are four 
branches (m = 1, . . . ,4 from bottom to top) with four 
states each. Also shown in Fig. || are the ipip* states for 
N = 64 (small circles). The solid lines are inferred from 
tpip* data for N — 2048 and represent the boundaries of 
the ijjtjj* continuum in the limit N — > oo. 

Why have we chosen to interpret the small green circle 
in Fig. |^ as an antipsinon and not as a magnon? Either 
choice is valid but we must heed the fact that antipsinons 
and magnons live in different physical vacua. 

When we interpret the small green circle as a magnon, 
then it coexists in the magnon vacuum with a macro- 
scopic number of fellow magnons (small black circles). 




FIG. 9: Energy versus wave number of all tptp* scattering 
states at M z /N = j and q > for N = 64 (o) in comparison 
with the corresponding free ipip* states (+). Also shown are 
the tptp* states for iV = 16 (large red circles). The inset shows 
the energy-momentum relations of the psinon (0 < \p\ < 7r/4) 
and the antipsinon (n/4 < \p\ < 37r/4). 



The collective excitations must then be viewed as con- 
taining a finite density of magnons (for N — > oo), in 
which the magnon interaction remains energetically sig- 
nificant for scattering states. The nonvanishing interac- 
tion energy obscures the role of individual magnons. 

On the other hand, when we interpret the small green 
circle as an antipsinon, then it lives in the psinon vacuum, 
i.e. almost in isolation. The only other particle present is 
a psinon (large green circle). In the limit N — > oo, the in- 
teraction energy in a psinon-antipsinon (ipip*) scattering 
state vanishes. 

The energy-momentum relations of the two quasi- 
particles (see inset to Fig. ^) can be accurately inferred 
from N — 2048 data for the spectral thresholds of the 
tj}\jj* states. The psinon dispersion e^{p) is confined to 
the interval at < \p\ < 7r/4 (solid line) and the antipsi- 
non dispersion e x p*{p) to 7r/4 < \p\ < 37r/4 (dashed line). 
The lower boundary of the iptp* continuum is defined by 
collective states in which one of the two particles has zero 
energy: the antipsinon for < \q\ < n/2 and the psinon 
for 7r/2 < \q\ < n. The upper boundary consists of three 
distinct segments. 

For < q < 0.3935 the highest iptp* state is made up 
of a zero-energy psinon with momentum = and 
an antipsinon with momentum — Tt/4 + q. Likewise, 
for 37r/4 < q < 7r, the states along the upper continuum 
boundary are made up of a maximum-energy antipsinon 
(with = 37r/4) and a psinon with p^ — — 37r/4 + q. 
In these intervals, the shape of the continuum boundary 



is that of e^j, or e 



tip- 



When the two delimiting curves are extended into the 
middle interval, 0.3935 < q < 37r/4, they join in a 
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cusp singularity at q = ir/2. Here the maximum of 
e ip{Pip) + e ip* iPii* ) subject to the constraint +p^* = q 
does not occur at the endpoint of any quasi-particle dis- 
persion curve. Consequently, the iptp* continuum is par- 
tially folded about the upper boundary along the middle 
segment as is evident in Fig. |. 

The (+) symbols in Fig. | represent N = 64 data of 
free i^nj}* superpositions generated from the ^ and tp* 
energy-momentum relations. The vertical displacement 
of any (o) from the associated (+) thus reflects the inter- 
action energy between the two quasi-particles in a tptp* 
state. A comparison of the N — 64 data in Figs, [j] and ^ 
shows that for the most part the tp^* interaction energy 
is smaller than the tptjj interaction energy (Problem 4b). 

The lower boundary of the V^A* continuum touches 
down to zero frequency at q = and q = q s = 7r/4, 
which is a special case of Eq. (|l2|). Between q s and ir, it 
rises monotonically and reaches the value E — Eq = h for 
N — ► oo (Problem 8). A direct observation of the zero- 
frequency mode at q s was mado-4n a neutron scattering 
experiment on copper benzoateJlj 

When we lower M z , the soft mode at q s moves to the 
right, the number of 2m-psinon branches contributing to 
the ipip* continuum shrinks but each branch gains addi- 
tional states. At M z = 1 we are left with one 2-psinon 
branch extending over the interior of the entire Brillouin 
zone. This branch is equal to the lowest branch of 2- 
spinon states with dispersion e^q), Eq. (||). At M z = 
the i/jtj}* excitations disappear altogether. 

When we increase M z toward the saturation value, the 
zero-frequency mode moves to the left, and the num- 
ber of 2m-psinon branches increases but each branch be- 
comes shorter. At M z = N/2 — 1, the two-parameter set 
collapses into a one-parameter set consisting of one 2m- 
psinon state each for m = 1,2,... , N/2 — 1. In part I 
these states were identified as 1-magnon excitations with 
dispersion ei(<z) = J(l — cosq) (Problem 9). 



VI. 



LINESHAPES 



What have we accomplished thus far and what remains 
to be done? We have identified the model system and the 
dynamical quantity which is relevant for the interpreta- 
tion of inelastic neutron scattering data (Sec. I). We have 
configured the ground state at h ^ as the vacuum for 
psinon quasi-particles and related it to the ground state 
at ft. = 0, the spinon vacuum (Sees. || and III). We 
have introduced a method of calculating matrix elements 
for the dynamic spin structure factor S zz (q,Lo) via the 
Bethe ansatz (Sec. IV). We have it to identify among 
all the 2m-psinon states one continuum of collective ex- 
citations which contributes most of the spectral weight 
to S zz (q,uj): the ^ip* states (Sec. [v]). Here we use the 
spectral information and the transition rates, all evalu- 
ated via Bethe ansatz, to calculate the lineshapes relevant 
for fixed-q scans in the neutron scattering experiment. 

To finish the task, we exploit key properties of transi- 



tion rates and densities of states of sets of excitations that 
form two-parameter continua in (q, w)-space for N — > oo. 
The two parameters are quantum numbers of the tp and 
ip* quasi-particles, e.g. the positions v, v* on the Jj- 
scale of Fig. || of the two green circles. For N — ► oo, 
the scaled quantum numbers v /N \v* /N become piece- 
wise smooth functions of the physical parameters q, u). 
The ipip* transition rates (scaled by N) then turn into 
a continuous function Mf/f (q,co) and the iptjj* density 
of states (scaled by TV -1 ) into a continuous function 
Dff(q,Lj). The ipip* spectral -weight distribution is the 
product Sff (q, lo) = Dff (q, cu)Mff (q, u) . 

In the following, we consider the case q = tt/2, where 
the tptp* continuum is gapless. The density of tpip* states 
is generated from N = 2048 data of 



Dff(q,u u ,) 



2ir/N 



(21) 



where v* — m marks the antipsinon quantum number 
in the tptp* continuum. The psinon quantum number 
v is adjusted to keep the wave number q of the ij)tp* 
state fixed. The resulting ordered sequence of levels sub- 
stituted into ( |2l| ) yields the graph shown in Fig. |lC|(a) 
(Problem 3b). 

The function Dff (tv/2,u> v *) rises from a nonzero 
value at u = slowly up to near the upper boundary, 
where it bends into a divergence. A log-log plot of the 
data near the upper band edge reveals that it is a square- 
root divergence. Finite- N data over a range of system 
sizes for the scaled transition rates 



Mtf(q,u„*) = N\(G\S*\v*)f 



(22) 



with q = tt/2 are shown in Fig. |To|(b). They confirm 
the smoothness of Mf/* (tt/2, ui). It has a monotonic 
w-dependence with a divergence at lo — and a cusp 
singularity at the upper boundary lujj ~ 1.679 J of the 
ijjip* continuum. 

The product of the transition rate function and the 
(interpolated) density of states is shown in Fig. |To|(c) . 
The curve fitted through the data points represents the 
ipip* lineshape of S zz (q,uj) at q = tt/2. Its most dis- 
tinctive feature is the two-peak structure caused by di- 
vergent transition rates and a divergent density of states 
at the lower and upper band edges, respectively. The 
divergence at lo = is a power law, Lu~ a , with an expp-. 
nent that is exactly known from field theoretic studies:E-3 
a = 0.4688 .... The strength of the divergence at 
uiu depends on whether or not the cusp singularity of 
Mff (it/2,uj) starts from zero at lo ~ lojj (Problem 10). 

The observability of this lineshape in a neutron scat- 
tering experiment also depends on the relative i/ji/j* con- 
tribution to the integrated intensity of S zz (tt/2, lo). With 
the tools developed here, we can show that the i/ji/j* states 
contribute at least 93% of the spectral weight at this 
particular wave number (Problem 11). A more complete 
analysis of the lineshapes and the integrated intensity 
of S zz (q,Lo) for the model system at M z /N — j can be 
found elsewhere. Ii3 
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FIG. 10: (a) Density of iptp* states at q = 7r/2 evaluated via 
(|l|) from Bethe ansatz data for N = 2048. (b) Transition 
rates ( |22[ ) between the psinon vacuum and the ipip* states 
at q = 7r/2 for N = 12,16,20,24,28,32. (c) Lineshape at 
q — 7r/2 of the ipijj* contribution to S zz (q,cj). All results 
pertain to M„/N = \. 



Where do we go from here with this series of tutorial 
papers on the Bethe ansatz? In the next installment, we 
plan on making the exchange interaction in ([!]) uniaxi- 
ally anisotropic. This well set the stage for the investiga- 
tion of the following topics. Like the magnetic field, the 
anisotropy is a useful continuous parameter which brings 
about interesting effects in the excitation spectrum of the 
spin chain. Unlike the former, the latter does affect the 
wave functions of the system. The Bethe ansatz offers 
us a close-up look into these changes. Foremost among 
them is the transformation of 2-spinon scattering states 
into 2-spinon bound states and the consequent changes 
of their roles in dynamic structure factors. 



VII. PROBLEMS FOR FURTHER STUDY 

(1) (a) Show that if you subject {Sf (t)S$(0)}, 
fi = x,y,z, at T = for a cyclic chain of N 
sites to a space-time Fourier transform, S^^{q,to) — 

Ed- e t,(i "° J_ + r ^e iut (^'(«)^(0)), you end up 
with expression (||). (b) Show that the static structure 
factor S l _ lfi (q) = (G\S£ S^ q \G) as obtained from a single 
diagonal matrix element is equal to the integrated inten- 
sity of @, f+™(du/2TT)S^( q ,Lo) = £„J(G|S£|m)| 2 as 
obtained from a sum of off-diagonal matrix elements. 

(2) Consider the lowest-lying 2-spinon excitation at 
q = 7r/2. Find the pattern of its Bethe quantum numbers 
Ii for arbitrary N. Use Eq. (119) to calculate the excita- 
tion energy E^^N) of that state for a range of system 
sizes. Plot the spinon interaction energy E 7r /2(N) — TrJ/2 



versus 1 /N to verify the TV-dependence suggested in the 
text. Repeat the same task for the lowest 2-spinon exci- 
tation at q — 7r r _IIere the reference energy is zero. Verify 
the exact resulO E„(N) ~ a/N, a — 7r 2 /2 for this exci- 
tation via extrapolation. 

(3) (a) Use the result (1130), e(q 7 q) — irJ\ sin ^qcosq\, 
< q < 9j f° r the 2-spinon spectrum at some wave num- 
ber < q < 7T to calculate the 2-spinon density of states 
analytically via D{q,uj) = J Q 9 dq8 (uj — e(q, q)). The re- 



sult D(q,ui) 



M) 



-1/2 



for CL{q) < ueu{q) is 



then to be multiplied by the exact transition rate function 



M(q, uj) derived in Ref. 10 to produce the exact 2-spinon 
part of S zz {quj). (b) Reproduce the analytic result for 
D(tt, uj) computationally from finite- N data via Eq. ( pi] ) 
with v* now labelling the 2-spinon states at q = n in 
order of increasing energy. 

(4) (a) For sufficiently large N the psinon interac- 
tion energy in the 2-psinon scattering states is of the 
form AE^(q) ~ e^(pi,p 2 )/N, where q = pi + p 2 
and e^,(pi,p2) depends smoothly on the psinon mo- 
menta p\,p2- Since we do not know the psinon energy- 
momentum relation e^(p) analytically, use the 2-psinon 
lower boundary for N* > iV as an approximation for 



e4p±n/4) in AE^\q)^ E^\q)- 



■e^(pi)-e^(p 2 )- Plot 

N/S.E^l\q) versus (pi,P2) for 2-psinon states of systems 
with N* 3> N 1 judiciously chosen. Compare the 
properties of the newly found function e^{px,P2) with 
those of the function e sp (p\,p2) established in Fig. || for 
spinons. (b) Carry out the same procedure for the tpip* 
states. The lowest branch of states is the same as in (a). 
Compare the trends in quasi-particle interaction energies 
with increasing energy of the 2-psinon states and ij)tp* 
states. 

(5) Identify a complete set of generators \j)o,j = 
1, . . . , d for the ( |l~5| ) in the subspace N = 6, r — 3. Verify 
that the dj associated with these generators add up to 
D. Establish the sets Jk for k = 0, tt/3, . . . , 5tt/3. Show 
that the add up to D. Calculate the Bethe coefficients 
dj of the d — 4 generators for all D = 20 solutions of the 
Bethe ansatz equations by using the data from Table IV 
of part I. Show that all dj with j (jt vanish. 

(6) Use the results of Problem 5 to calculate (for 
N = 6, r = 3) the structure factor (G\S%S» q \G) and 
the transition rates \(G\S^\m)\ 2 . Show that these data 
satisfy the general result of Problem 1 (b) . 

(7) Show that \{G\S* =0 \G)\ 2 = M|/iV and use this 
result to test your computer program which calculates 
transition rates from Bethe wave functions. 

(8) The 2m-psinon branch (0 < m < M z ) of the ipip* 
continuum at q > is specified by the r = N/2 — M z 
Bethe quantum numbers I\ — — N/4 + M z /2 + | — m, 
-N/4 + M z /2 + \ < I 2 < . . . < I R < N/4 - M z /2 - \. 
(a) Identify the Bethe quantum numbers of the ipip* state 
with the lowest energy above the ground state \G). Show 
that the wave number of this state is k — kc = 7r(l — 
2M Z /N). Show by numerical extrapolation of finite- N 
data that its excitation energy tends to zero as N — > oo. 
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(b) Identify the Bethe quantum numbers of the sole ipip* 
state with q = ir — 2-k/N . Plotting the excitation energy 
of this state versus M z /N yields a set of data points that 
converges, as N — > oo, toward the inverse magnetization 
curve h(M z ). 

(9) Use the exact 1-magnon wave functions (16) with 
Bethe coefficients (18) in expressions @-@ to calcu- 
late the iptp* transition rates |(G|S||A)pTbr M z = N/2- 
1 and arbitrary N. Note the different TV-dependence for 
q = and q ^ 0. Show that the relative tptp* spectral 
weight in S zz (q,u>) is 100%, but the absolute intensity 
for q ^ is only of 0(N- r ). 

(10) (a) Use the lowest-energy excitation for N = 
12, 16, . . . , 32 to perform a nonlinear fit of the expression 
Mj/ (ir /2,o;p) ~ ai + a20J~ a and compare the exponent 
value thus obtained with the field-theoretic prediction 
quoted in the text, (b) Use all data at huj/J > 1 to fit 
the expression Mf^(n/2,w p ) ~ h + b 2 {uJu - u) . Per- 
form an alternative fit in which b\ is forcibly set equal to 
zero. 

(11) For fixed N, the integrated intensity is the expec- 



tation value (G\S*/ 2 S*/ 2 \G) and the ipip* part thereof the 



sum of transition rates Ylu'=i I (^1^/21^*) I 2 as wor ked 
out in Problem 1. Evaluate these quantities from the 
Bethe wave functions for N = 12, 16, . . . and extrapo- 
late the ratio find the relative spectral weight of the tpip* 
excitations in S zz (ir/2,u>). 
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